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ABSTRACT 



Aims. We have examined a stream-source model for the production of the cometary tails observed in the Helix Nebula NGC 7293 in which a 
transonic or moderately supersonic stream of ionized gas overruns a source of ionized gas. We have compared the velocity structures calculated 
with the available observational data. We have also investigated the suggestion that faint striations visible in the nebular gas are the decaying 
tails of now destroyed cometary globules. 

Methods. We have selected relevant results from extensive hydrodynamic calculations made with the COBRA code. 

Results. The velocities calculated are in good agreement with the observational data on tail velocities and are consistent with observations 
of the nebular structure. The results also are indicative of a stellar atmosphere origin for the cometary globules. Tail remnants persist for 
timescales long enough for their identification with striations to be plausible. 
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1. Introduction 

Many PNe contain extended cometar y tail-like structures. 
Thes e emanate from dense, neutral (Meaburn et al 1 119921: 
iHuggins et all Il992l) globules with radiatively ionized sur- 
faces and point away, comet-like, from the central star. The 
most studied are those in the highly evolved PNe NGC 7293 
(the Helix Nebula). A brief review of the two broad classes 
of models proposed to account for them is given by Dyson 
(2003). "Shadow" models (e.g.. Ivan Blerkom & Arnvlll972t 
ICanto etafl Il998( lLopez-Martin et al.ll200ll) use an opaque 
cloud to shield gas from the direct radiation field of the PNe nu- 
clear star. The shadow region is energized by a diffuse radiation 
field of different quality from t he dire ct stellar field. "Stream- 
source" models JDvson et alJ i 19931) - henceforth Paper I, 
iFalle etalJd2002l) - henceforth Paper II) are based on mass in- 
jection from embedde d clumps in global flows. T here are three 
relevant length scales jHartauist & Dvsonll993l) . The smallest 
is the mass injection scale and the largest scale is that for as- 
similation of injected gas into the global flow. The intermediate 
scale where injected material largely retains its own identity, 
produces a tail feature that is dynamically affected by momen- 
tum transfer from the global flow. 

iLefloch & Lazarefll(ll994 argued for tail formation during 
the photoionization of dense globules, a suggestion also made 



Send offprint requests to: J.E. Dyson 



by iRiikhorst et alJ J2006I) who did not, howev er, present any 
models. Yet as shown bv lPavlakisetalJfeOOlll. the diffuse ra- 
diation field which was neg lectedbv lLefloch & Lazarel ( ll994 
can suppress tail formation if it has an intensity of a few percent 
of the direct stellar flux. Given that there are several thousand 
tails distributed throughout the nebula with concommitant vari- 
ations in direct flux, diffuse flux and inevitably globule proper- 
ties, it is hard to see how the tail system could be produced by 
this mechanism. Stream-source models also have the advantage 
that they are inherently dynamic in distinction to the majority 
of shadow models where the dynamics of the g as is neglected, 
though the model of lLopez-Martin et alJ J200ll) includes a flow 
of ionized gas perpendicular to a tail feature as the gas exits 
from a D-critical ionization front at the ionized sound speed. 
However, the only detailed o bservational study of th e velocity 
structure of the Helix tails by iMeaburn et alJ Jl998t) shows no 
evidence for significant ionized gas velocities perpendicular to 
the tails. 

The ratio of tail length to head diameter of the Helix 
cometar y features ranges from w 1 to * 100 in the spectacular 
knot 38 jMeaburn et alj[l998l) which has a deprojected actual 
tail length of ta 2.5 x 10 17 cm. Tails are aligned to within a 
few degrees along radial vectors originating at the central star. 
Our discussions are directed specifically at the Helix tails and 
clumps and their generalisation to other PNe may be inappro- 
priate. 



2 



J.E. Dyson, J.M. Pittard, J. Meaburn and S.A.E.G. Falle: The tails in the Helix Nebula NGC 7293 



In this paper we re-examine the stream-source models pro- 
posed in Paper I using extensive calculations of stream-source 
interactions by Dyson et al. (in preparation - henceforth Paper 
IV). In Paper I, a hot subsonic stream of bottled-up shocked 
stellar wind gas was assumed to interact with a source of cooler 
gas to produce tails. Serio us object i ons to this tail generation 
mech anism were raised bv lO'Delll lEo 00) and Meaburn et al. 
J2005I) . The main problem is that the central cavity in the Helix 
is occupied by a hot (T e > 2 x 10 4 K) He + emitting gas and it 
is hard to see whe re the hot shocked stellar w ind gas resides. In 
fact much earlier. IMeaburn & White! Jl982l) had raised a simi- 
lar objection because of the presence of an inner [OIII] emitting 
shell that would shield the cometary globules from the direct 
impact of any fast (subsonic or supersonic) wind. 

Therefore we examine the suggestion that the cometary 
tail heads are over run by a transonic or supersonic stream of 
photoionized gas JMeaburn et alJ fl992l Il998t lO'DelllkOOol 
IMeaburn et al J 12005 ). We compare result s with the knot dy- 
namics observed by Meaburn et alJ Jl998l) and show these are 
reasonably well reproduced. We briefly discuss the significance 
of clump motions for the se models. 

IMeaburn et alJ J2005I) noted the existence of "headless" ra- 
dial spokes close to the central star of NGC 7293 and suggested 
that these spokes were tails left be hind after the i r com etary 
knots were photoe vaporated away. IHenrv et"afl J 19991) and 
lO'Dell et alJ j2004 also noted similar striations in the nebular 
gas. We examine the effects of switching off the mass injection 
and show that tails persist for typical timescales of thousands of 
years before they are advected away by the stream flow, le nding 
plausibility to the suggested origin JMeaburn et all2005l) . 

2. Stream-source tail production 

Paper I argued that long thin tails are produced only when sub- 
sonic streams interact with subsonic mass injection (all Mach 
numbers are defined in the reference frame of the clump). 
Paper II dealt with the interact ions of hypersonic an d transonic 
streams with injected matter. IPittard et alJ J2005I henceforth 
Paper III) investigated the interaction of hypersonic flows with 
multiple clumps and the interaction of transonic flows with two 
adjacent clumps. We will consider only single clump interac- 
tions although both photo graphs of tails and the recent estimate 
of 23000 cometary knots JMeixner et al.ll2005T) suggest multi- 
ple clump interactions might occur in the Helix. 

In Papers I to III it was assumed that there was a large tem- 
perature contrast between a hot global flow and injected gas 
with the stream gas behaving adiabatically and the injected gas 
behaving isothermally. Tail dynamics and morphology are de- 
termined by the thermal behaviour and temperature contrast of 
the stream and source gas as well as by the Mach number of the 
incident stream. Paper IV investigates a broader range of inci- 
dent Mach numbers and temperature contrasts using spherical 
clump geometry (Paper III shows that the differences between 
spherical and cylindrical geometry are small). Results are also 
given for the case where both stream and source gases behave 
isothermally. This is appropriate if the stream gas as well as 
the source gas is photoionized and we utilise these results here. 
Computational details are given in Papers III and IV. 



An important result from these new calculations is that 
identifiable long thin tails persist to appreciably higher Mach 
number interactions than suggested in Papers I and II. This is 
quantified in Paper IV. 

Mass loss from the clumps in the Helix is due to photoion- 
ization. Clump images show the mass injection is strongest in 
the direction towards the central star. In Paper III it was shown 
that anisotropic mass injection made very little difference to the 
flow morphology in the case of a hypersonic wind interacting 
with source material. Comparison of tail structures for subsonic 
and transonic stream interactions where there is anisotropic in- 
jection is given in Paper IV. Results for a mass injection rate 
contrast of 20: 1 (front to rear) show small differences in the tail 
morphologies and tail velocity structures between the isotropic 
and anisotropic cases once the flow has gone a few injection 
radii. We use results for isotropic mass injection here. 

Although the models of Paper IV are far more extensive 
than previously given, they have still a major simplification. 
They contain only two gas phases (stream and injected gas) 
which may mix. The Helix tails possibly represent a three 
phase system that includes io nised stream gas and both neutral 
and io nised injected gas fe.g.. lO'Dell & HandroiJfggcllO'Deri 
l2000l) . To treat this properly requires the inclusion of radiation 
transfer. 

3. Tail production 

IMeaburn et ail Jl998l) . lO'Delll J2000I) and IMeaburn et ail 
J2004 have shown that it is likely that clumps are overrun by 
photoionized gas that may be ionized AGB wind possibly con- 
taminated with evaporated knot gas. We consider tail formation 
in an interaction where photoionized gas overruns clumps that 
lose material by photoionization. An isothermal-isothermal as- 
sumption for the thermal behaviour of the stream and source 
gas is appropriate. 

It is unlikely in general that the clumps are being over- 
run by the very hot He ++ gas since the clumps have an aver- 
age outward velocity of 14 kms 1 JMeaburn et alJI 19981) and 
the hot gas is dynam ically inert with an e xpansion velocity of 
less than 1 1 km s _1 JMeaburn et all2005l) . On the other hand, 
there might be clumps that are overrun by this hot gas since 
the velocity dispersio n of the knots is around 6-8 kms 1 
JMeaburn et aljfl998l) . The temperature c ontrast betwee n the 
stream gas and injected gas may be high. lO'Delll J1 9981) sug- 
gests that the temperature of the He ++ gas is i n excess of 
2 x 10 4 K (because of the lack of ++ cooling'). IHenrv et alJ 
Jl999l) have modelled the ionization structure of the nebula and 
shown that the gas temperature drops from about 4x 10 4 K near 
the star to 10 4 K at a distance where helium is about 50% in 
the form of He ++ . 

IMeaburn et alJJl998l) and lMeaburn et al.lJ2005l) favour the 
overrunning of clumps by expanding photoionized gas in the 
region where IOIII 1 is emitted . This gas has a temperature of 
about 10 4 K JHenrv et al.ll999t) . and is overrunning the clumps 
with a velocity of about 17 km s _1 relative to the clump global 
expansion velocity of 14 kms"' . 

The temperature of gas injected from the source is uncer- 
tain. [o^Del^uiO H(X)il) note that the temperature of photoion- 
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Fig. 1. Tail density structures for sound speed ratio cs/q =2:1. 
The stream Mach number, M = 1, 2, 4 for the left, middle, and 
right panels respectively. Each panel shows a region < r < 
30, -100 < z < 10. The logarithm of mass density is shown, 
with a scale of -0.69 to 0.39 (left), -0.97 to 0.99 (middle), and 
-0.99 to 1 .60 (right). Blue denotes the lowest density, while red 
denotes the highest. The stream density is 0.25. An injection 
region of unit radius is at the origin. See Paper II for further 
details about the computations. 



ized gas flowing from a clump towards the star decreases out- 
wards from the density peak just behind the ionisation front 
on the clump surface. Some injection will occur from the sides 
of clumps where the diffuse radiation field softens the average 
photon energy leading to gas temperatures that could be a fac- 
tor of around 2 lower than the temperature of gas excited by the 
direct radiation field ( I Van Blerkom & Arnvll 19721 ICanto etafl 
1998). 

As a reasonable range for the temperatures we use calcu- 
lations from Paper IV for ratios of isothermal sound speeds 
(cs/q, stream gas to injected gas) of 1:1 (equal temperatures) 
and 2:1 (a temperature ratio of 4:1). In Fig.^ we gi ve tail den- 
sity structures for a sound speed ratio of 2: 1 for incident stream 
Mach numbers M = 1,2,4. 

In Fig. |3 we give the same information for a sound speed 
ratio 1:1. However a sound speed ratio of unity gives a tail den- 
sity effectively the same as that of the stream and the plots of 
Fig. |3 therefore do not show systematic tail structure. The re- 
sults of Paper IV demonstrate that the injected gas maintains 
its identity and there is very little mixing between it and stream 
gas. So to show this, we plot in Fig.[3]the advected scalar that 
distinguishes between the stream and injected gases. A recog- 



Fig. 2. As Fig.^but for equal sound speeds with M = 1,2,4 
(left, middle, and right, respectively). The density scale extends 
from -0.52 to 0.39 (left), -0.97 to 0.99 (middle), and -1.02 to 
1.59 (right). The stream density is 1.0. 




Fig. 3. As Fig.|2]but showing advected scalar structures. Blue 
represents the stream gas, while red represents injected material 
(i.e. material evaporated from the clump into the flow). 
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nisable tail of injected gas is embedded in stream gas. Some 
difference in emission or absorption properties (e.g., a different 
gas-to-dust ratio) is needed to make the tail visible as is clear 
from the density plot but it is a dynamical possibility for tail 
production. 

4. Tail velocity structures 

Important diagnostic informat ion comes from the tai l velocity 
structure in knot 38 given by iMeaburn et al.l 0998I) . In view 
of our model limitations, we use mass-weighted velocities for 
comparison with the data. 

In Fig.|4]we show mass-weighted velocity profiles along 
tails for several incident Mach numbers and two stream-source 
gas temperature ratios. We give results only for the transonic 
and supersonic isothermal-isothermal flows that are relevant to 
the present paper. 

For streams that are roughly in excess of transonic, the ve- 
locity increases quickly over a few injection radii. The initial 
increase is followed by a velocity decrease for Mach numbers 
less than about 3 then the velocity rises to a roughly coasting 
velocity. Once the Mach number is around 3, the initial increase 
is just followed by a very gentle rise or even roughly constant 
velocity over distances of tens of injection radii. Appreciably 
greater distances are probably not relevant to the Helix tails 
where the very longest tail (knot 38) has a length about 100 
times its head diameter and we display results accordingly. 

The final velocity reached depends on the temperature con- 
trast between the stream and source gas and is higher (in terms 
of the stream velocity, Vs) the lower the temperature contrast. 
This is as expected since the temperature contrast determines 
crudely the stream-injected gas density contrast and there is 
more momentum transfer from unit mass of the stream to unit 
mass of injected gas the less this contrast. At sufficiently high 
incident Mach number and low enough temperature contrast, 
the tail material achieves supersonic velocities in terms of the 
sound speed in the stream and obviously as well that of the 
injected gas (since the injected gas is either at the same tem- 
perature or cooler than the stream gas). 

To compare the tail models with the observational data we 
assume that a stream of velocity Vs interacts with a clump 
which has a velocity Vc, where both velocities are measured 
in the stellar frame. The relativ e velocity is therefore (Vs - Vc). 
We concentrate on knot 38 of IMeaburn etall Jl99«h . The gas 
at the en d of this tail has a velocity V r = 22 kms 1 relative to 
the head ( IMeaburn & Redmanl2003l) . Since this knot has a tail 
length-to-head diameter ratio of around 100:1 we take this to 
be equal to the coasting velocity from our calculations reached 
at 100 injection radii. 

We can write V T = P(M)(V S - V c ) = f3(M)Mc s where 
f3(M) is a constant dependent on the Mach number M of the 
interaction and the stream-source sound speed ratio and M 
( Vs - Vc)/cs. In Table^we give values of f3(M) for the cases 
considered here. 

In Table |2] we give tail velocities calculated for various 
Mach numbers for the two temperature ratios above using the 
global knot velocity of 14 kms 1 . We use a sound speed of 
11.7 kms~' (i.e. c; = 10 4 K) for the source gas. We also list 



Spherical clump, M = 1 .0, Isothermal-Isothermal 
c/c 5 = 1 .0 



Spherical clump, M = 2.0, Isothermal-Isothermal 



!(rj 

Spherical clump, M = 4.0, Isothermal-Isothermal 





c/c s = 1 .0 






c/c s = 0.5 











Fig. 4. Vz vs z for M — 1,2,4, and sound speed ratios 1 and 2 
(i.e. temperature ratios 1 and 4). The injection radius is r c . 

Table 1. The value of /3(M) as a function of M and cs/cj. 



M 


Cs/Cj = 1 


cs/ci = 2 


1 


0.91 


0.76 


2 


0.81 


0.54 


4 


0.66 


0.43 



values of Vs, the stream velocity in the stellar frame, for com- 
parison with the value of 31 km s _1 of the overrunning stream 
found bv lMeaburn et al.Nl998l) . 

The observed acceleration towards positive radial velocities 
along the tail from kn ot 38 is shown particula rly well in the line 
profiles in Fig. 13 of IMeaburn et alJ Jl998l) and in the veloc- 
ity images in Figs. 4 and 12 of that paper (though first correct 
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Table 2. Tail ( V T ) and stream ( Vs) velocities ( km s 1 ). 





Cs/Ci 


= 1 


Cs/Cj 


= 2 


M 


K T 


V s 


V T 




1 


11 


26 


18 


38 


2 


19 


37 


25 


61 


4 


31 


61 


40 


107 



the error in the legend for Fig. 12 which includes four panels 
whose indicated positions should all be reversed i.e. bottom- 
left, bottom-right, top-left, top-right should read top-right, top- 
left, bottom-right, bottom-left respectively). With this correc- 
tion it can be seen in Fig. 12 that the tail furthest from the knot 
becomes more prominent in the more positive radial velocity 
range. 

We see that the tail and overrunning stream velocities are 
reasonably well reproduced with M « 2, cs/c; = 1 and M « 1, 
c s /cj = 2. 

Tails should show more complex velocity structure closer 
to the injection zone particularly with anisotropic mass injec- 
tion but we do not consider this further here. 

The setting up time for the tail of knot 38 ~ 2 - 4 x 10 3 yrs 
with a length of 2.5 x 10 17 cm and a tail gas velocity of 
10 - 20 kms _I . This i mplies formation after fast win d activ- 
ity has ceased. Although lPatriarchi & Perinottollll99lh showed 
from IUE observati ons that 60% of central s t ars of PNe emit 
fast particle winds, ICerruti-Sola & Perinottol dl985l) failed to 
detect one from the central star of NGC 7293. Shorter tails 
take proportionally shorter times to form and more "typical" 
tails with lengths of a few times 10 15 cm could be set up in 
timescales of the order of a hundred years. It is not clear why 
the tails are generally fairly short. One possibility might be time 
fluctuations in stream or source injection properties. 

5. Clump motions 

The effects on tail production mechanisms of systematic clump 
motions, radial or non-radial, have received little attention. The 
non-radial component of the Helix knot expansion velocity is 
unclear. On general grounds, non-radial velocity components 
of clumps seem inevitable if clumps form in dynamical in- 
stabilities. If they are generated by shell fragmentation, they 
should have some dispersion around the shell velocity at frag- 
mentation since it is highly unlikely that shells are strictly 
sp herical. 

Diamond & Kemballl d2003l) measured global proper mo- 
tions of the SiO maser spots in the Mira variable TX Cam 
that correspond to velocities of about 10 kms 1 and noted 
that "significant" non-radial velocity components are present. 
These motions are in a stellar atmosphere of radius ~ few x 
10 13 cm. If these spots move out to radial distances ~ 10 17 cm, 
conservation of angular momentum ensures very small non- 
radial velocity components («: 1 kms~' ). 

In Paper IV results are given for situations where clumps 
move at an angle to a global flow. To a good approxima- 
tion, well-developed tails lie along the resultant velocity vec- 
tor of the stream-clump velocities. If we assume that the 



Table 3. V tan ( km s 1 ) for an alignment angle of 5°. 



M 


Cs/Cj = 1 


cs/t'i = 2 


1 


1 


2 


2 


2 


4 


3 


4 
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stream is strictly radial, then tails are aligned to an angle 
a a arc tan[V tan /(Vs - Vc)] with the radial direction, where 
V tan is the non-radial velocity component of the clump veloc- 
ity. In Table [3] we give maximum values of V tan that allow an 
alignment angle of 5° to the radial direction where we take this 
angle as characteristic. 

The best-fit cases for Vt and Vs detailed above require very 
small values of V tan (< 2 kms" 1 ) for alignment. This is un- 
likely to be satisfied with shell instability formation but is more 
consistent with the much lower non-radial velocity components 
expected if the clumps originated in stellar atmospheres. 

6. Tail persistence 

We have examined the persistence of tail s where the injection is 
switched off to quantify the proposal bv lMeaburn et alJ d2005l) 
that this is responsible for the production of headless spokes. 
In Fig. |5] we show time frames for the evolution of the tails 
produced for the two best fit cases in Section |4] (i.e. M ta 2, 
cs/ci = 1 and M as 1, cs/cj = 2). As before, for the case of 
a sound speed ratio of unity we give both the advected scalar 
plots and the density plots. Tails persist over timescales of at 
least 50 or so (in dimensionless units). The actual time unit 
is equal to 2rc/VF, where Vf is the gas velocity in the clump 
frame (i.e. Vs - Vc). If we take the sound speed in the injected 
gas to be 11.7 km s~ 1 as before, the flow speed in both cases is 
23.4 kms -1 . If we take the injection radius as about 10 15 cm 
corresponding to a typical clump diameter, the tails persist for 
times at least ~ 1400 yr. It appears plausible therefore that 
headless tails are the remains of tails behind recently destroyed 
clumps. 

7. Conclusions 

Linear structures have been observed in several planetary neb- 
ulae. This morphology has led to two classes of models. 
However morphology itself is not sufficient to discriminate be- 
tween them. Velocity structure is the key and only for two knots 
(14 and 38) in the Helix nebula is the necessary data available. 
The acceleration along the tail associated with the spectacu- 
lar knot 38 suggests strongly that there is momentum transfer 
between a directed flow and a mass source. The advanced evo- 
lutionary state of the central star, and the details of the nebular 
structure, rule out any hypersonic stellar wind-mass source in- 
teraction. Although attractive in many respects, the interaction 
with a subsonic hot stream of shocked stellar wind gas is also 
ruled out because the existence of a reservoir of hot gas cannot 
be reconciled with the nebular structure. Much more plausible 
is tail formation taking place when a transonic or moderately 
supersonic stream of photoionized gas overruns gas injected by 
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Fig. 5. Evolution of the M = 2, c,/c, = 1.0 isothermal-isothermal simulation after the source injection is switched off. t=0 (left), 
6.0, 13.4, 28.3, and 58.0 (right). The isothermal sound speed=1.0 for both the injected and ambient flow, and the ambient flow 
speed=2.0. The top panels show density plots, while the bottom panels show the tracer for injected material. The colour scaling 
of each density panel is different in order to best display the structure. From left to right, logio p = -0.97 to 0.99, -0.14 to 0.28, 
-0.38 to 0.15, -0.38 to 0.12, and -0.28 to 0.08. The density enhancement relative to the stream density (logio p = 0.0) declines 
as time evolves. 



the photoionization of cometary globules. The injected gas is 
accelerated by momentum transfer from the stream and the ve- 
locities that are generated are in good agreement between those 
measured for the likely overrunning stream and the tail gas of 
knot 38 in the Helix. Thus the general scheme for the formation 
of cometary globules involves flows (or winds) overrunning a 



slowly expanding system of pre-existing dense neutral glob- 
ules. At least for NGC 7293, the knotty nature of the ambient 
neutral shells associated with this object surely precludes any 
other possibility (instabilities in a global expanding shell for 
example). In fact the very strong radial alignment of the tails 
suggests that when overrun, the globules have a velocity that 
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Fig. 6. Evolution of the M = 1, cs/c; = 2 isothermal-isothermal simulation after the source injection is switched off. t — (left), 
10.9, 21.7, 43.4, 76.4, and 120.4 (right). The isothermal sound speeds 1.0 and 2.0 for the injected and ambient flow (respectively), 
and the ambient flow speed=2.0. The colour scaling of each panel is different in order to best display the structure. From left 
to right, log io p = -0.69 to 0.39, -0.70 to 0.11, -0.81 to 0.04, -0.81 to 0.03, -0.74 to 0.04, and -0.74 to -0.32. The density 
enhancement relative to the stream density (logio p = -0.60) declines as time evolves. 



can at most have a very small non-radial component. This is 
consistent with density structures that originate at radii far less 
that the present scale size of the nebula. Finally, it also appears 
likely that the widely distributed striations seen in the nebular 
gas represent the remains of tail structures generated by mass 
injection from now destroyed globules. 

Acknowledgements. JMP gratefully acknowledges support from the 
Royal Society. This research has made use of NASA's Astrophysics 
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